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Abstract 

We report on the observation of vortex formation in a Bose-Einstein condensate of ^'^ Kb atoms. 
Vortices are generated by superimposing an oscillating excitation to the trapping potential in- 
troduced by an external magnetic field. For small amplitudes of the external excitation field we 
observe a bending of the cloud axis. Increasing the amplitude we observe formation of a grow- 
ing number of vortices in the sample. Shot-to-shot variations in both vortex number and position 
within the condensed cloud are observed, probably due to the intrinsic vortex nucleation dynamics. 
We discuss the possible formation of vortices and anti-vortices in the sample as well as possible 
mechanisms for vortex nucleation. 
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I. INTRODUCTION 



Superfluidity is a remarkable signature of the quantum nature of a given physical system. 
This phenomenon is consequence of the macroscopic occupation of a single quantum state 
and reveals itself as a frictionless flow without dissipation. The superfluidity was discovered 
in 1938 by Kapitza [l| and, independently, by Allen and Misener in a sample of liquid 
^ife at a temperature of 2.17K. In particular, formation of quantized vortices is one of the 
most outstanding features of superfluids, and a considerable amount of work has been done 
on systems like ^He or mixtures of ^He/^He 3|. 

The achievement of quantum degeneracy in trapped dilute atomic gases has opened 
new directions in the study of superfluidity and vorticity. In fact, these quantum fluids are 
natural testing grounds for studying properties and phenomena related to superfluidity due 
to the unique possibility of controlling many of the parameters of the system in an easy 
way. In particular, the observation of quantized vortices in a BEC sl and, more recently, 
in a quantum degenerate Fermi gas jo], has been considered as an undeniable proof of their 
superfluid nature. 

Besides the big tunability of the different parameters of these systems, there are many 
advantages in realizing these studies in quantum degenerate gases. The possibility of ex- 
panding the quantum atomic fluid, together with the fact that the low density regime yields a 
quite large healing length which characterizes the vortex core size, make possible the direct 
observation of the vortex characteristic distribution, geometry and location using optical 
absorption techniques. 

The first experimental production of a vortex in an atomic BEC was performed by 
Matthews et al. js] using a phase engineering technique. In that case, rotation was dy- 
namically induced by driving a transition between two different hyperfine spin states of the 
trapped atoms, using a rotating field. Vortices have also been produced by focalization 
of a laser beam on the condensate moving faster than the critical velocity Stirring 
the condensate with a laser beam has allowed many different experiments involving vortex 
formation Q|, determination of the transferred angular momentum to the cloud and deter- 
mination of threshold-frequency during vortex excitation j^l. Rotating the condensate with 



an asymmetric trapping potential 



10| was the first demonstration of a purely magnetic exci- 



tation scheme quite analogous to the rotating bucket experiment with liquid He Sj. Also, a 
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dynamical instability ll| lead to the nucleation of vortices, with subsequent crystallization 



in a lattice configuration. Recently, vortex formation has been also observed by merging 
multiple trapped BECs [3] and by coherent transfer of orbital angular momentum from 



optical fields to the condensed sample 

n 



13| . Vortices in superfiuids continue to generate 
interest both theoretically [l^ and experimentally [l^. 

In this communication we show experimentally that vortices can be nucleated in a BEC 
when an oscillating spherical quadrupole field is applied to it. We are able to transfer several 
quanta of angular momentum to the condensate, evidenced by the appearance of many 
vortices in the sample. This technique is similar to the theoretical proposal by Mottonen et 

n n 

al. [16| and recently by Z.F. Xu et a/.p/7| to pump angular momentum in a BEC using an 
external magnetic field to imprint local Berry phase in the sample, though we do not believe 
that this is the actual mechanism in our case, mainly because in both protocols an inversion 



of the magnetic bias field is needed and that is not the case in our experiment. 

II. EXPERIMENTAL SETUP 

The experimental setup to produce the quantum sample is described in greater detail 



in Ref. 18|. Briefiy, we collect 10^ ^''Rb atoms at lOOfiK in a MOT. The atoms are then 
transferred in the |2,2) hyperfine state and then loaded into a QUIC-type magnetic trap 
[l^. Once atoms are in the QUIC potential, RF-induced evaporation is used to obtain 
quantum degeneracy. The trapping potential is harmonic with frequencies given by = ujq 
and Uy = Uz = 9u!q, where cjq ~ 27r x 23Hz. We extract the parameters of the atomic cloud 
(number of atoms, density profile and temperature) by imaging it on a CCD camera after 
free fall of 15 ms once the trapping potential is turned off. The imaging beam propagates in 
the xz plane, making a 45° angle with the loffe coil shown in Fig|T](a). We typically 

produce BEC samples containing 1 — 2 x 10^ atoms. 

After reaching BEC, an extra field, produced by a pair of anti-Helmholtz coils is super- 
imposed on the QUIC trap field as shown in Fig{T](b). The axis of the extra coils is close to 
the trap axis (the angle between these axes is 5 ~ 5°). The center of the extra field, defined 
by the zero-field amplitude position, is close to the QUIC trap minimum. An oscillatory 
current is applied to the coils, always having the same sign and always starting from zero, 
so we do not give an abrupt kick to the condensate in the beginning of the excitation phase 
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FIG. 1: Schematics of (a) above view, and (b) side view of the imaging beam and the coil arrange- 
ment in our system, showing the QUIC trap coils and the ACC set of coils. The drawing is not to 
scale and the ACC misalignment is made bigger for clarity. 



{Icoii = -^0 [1 — cos(f2t)]). The experimental sequence is then as follows: finished the evapo- 
rative cooling, the external field is turned on for tens of milliseconds. Once the external field 
is turned off, atoms are held in the static magnetic trap for extra 20ms prior to be released 
from the QUIC trap. Atoms are imaged after a 15ms time-of- flight. For this experiment we 
restrict ourselves to excitation periods running from zero to 60ms and maximum amplitudes 
of the spatial gradient of the field up to ~ 200 mG/cm. 

We indicate the set of excitation coils as AC-coils (ACC). The effect of field added by the 
ACC over the static trap is to produce a combination of translation, i. e. a displacement 
of the minimum trap, rotation, due to the angular dislocation of the coils, and a shape 
deformation, due to the asymmetry caused by the combination of the ACC field with the 
QUIC field. A simulation of the field for experimental parameters close to the real ones is 
shown in FiglS 
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FIG. 2: Simulated movement/deformation of the trapping potential under the influence of the 
external oscillatory field, (a), (b), (c) are the XY plane and (d), (e), (f) represent the XZ plane. 
Images (a) and (d) are at t=0, the starting point of the excitation, (b) and (e) are at midway 
between zero and the maximum amplitude of the external field {t = 7r/2a;) and (c) and (f) are at 
maximum amplitude field (t = tt/uj). The field of view is the same for all images. 

III. VORTEX OBSERVATION 

The spatial gradient of the field produced by the ACC is small, on the order of tens to 
hundreds of miliGauss/cm for the axial gradient, which in amplitude is comparable with the 
field variation experienced by the atoms in the BEG. We will show that for certain values of 
amplitude and frequency, this configuration transfers angular momentum to the BEG and 
hence produces vortices. 

A small amplitude of excitation (< 40 mG/cm) produces only a bending oscillation of 
the trapped cloud axis as shown in Fig|3l This tilting on the superfluid axis shows that 
excitation by the oscillatory field is able to mechanically transfer angular momentum to 
the atomic cloud. This effect can be related to scissor modes extensively investigated in 
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FIG. 3: Observation of cloud axis bending oscillations after low amplitude (< 40 mG/cm) applica- 
tion of the ACC field. Each image is the result of a different run of the experiment under similar 
experimental conditions and shows the optical density of the sample after 15ms of free expansion 




FIG. 4: Optical density of the sample showing a typical single vortex pattern after application of 
the ACC oscillatory field after 15ms of free expansion. 



references 



20 



2l|. 



The oscillatory excitation nucleates vortices only after reaching certain amplitude for the 
ACC field. For a fixed time of excitation (20ms), f2 = 27r x 200Hz, and amplitudes up to 
40 mC/cm of the axial field of the ACC coils, we observe no evidence of vortex formation. 
Nevertheless, the stirring of the whole cloud as in Figj3]is present. Increasing the amplitude 
further, the appearance of a single vortex type pattern starts to take place, as shown in FigJH 
Up to about 90 mC / cm we observe either zero or one vortex when the same conditions are 
considered. The successive repetition of the experiment under the interval between 40 and 
90 mC/cm reveals an average number of vortex formation of 0.6 ± 0.3. We call this as 
a single vortex zone of amplitude. For amplitudes larger than 90 mC/cm and up to 160 
mG/cm, most of the images show the presence of two or more vortex-type profiles, with a 
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FIG. 5: Typical atomic optical density images showing the appearance of (a) two and (b), (c) 
multiple vortices patterns after different excitation conditions. Images were taken after 15ms of 
free expansion. 



predominance of three vortices and occasional formation of one or more than three vortices. 
Typical images acquired in this regime are shown in FigJS] Again, several repetitions of 
the experiment within this amplitude range under equal experimental conditions allows to 
observe a multiple vortex zone, with an average number of vortices of 2.6±1.2. A considerable 
increase in the number of vortices is observed with larger oscillation amplitudes. We will 
not address this region in this communication since the vortices patterns observed are not 
regular nor the cloud seems to behave in the same way. These evidences might indicate that 
we access other regimes of vortex formation and equilibrium, not yet understood and out of 
the scope of this work. 

The variation of vortex number observed in our system when operating at equal conditions 
is not clear. It could be related to the mechanisms involved in the vortex formation and 
eventual decay during our excitation or it can be due to some out of control instability 
generated by the presence of the ACC field and its mutual inductance with the QUIC trap 
coils. Another important feature that has to be taken into account for explaining this 
fluctuation is that there can be a kind of energy degeneracy among different vortex states, 
resulting in equal probability of formation at equivalent experimental conditions. 

Figj6] summarizes our experimental observations and clearly shows that we increase the 
angular momentum transfer to the condensate as the amplitude of the excitation is increased. 
This is evidenced by a fast growing average number of vortices as the amplitude of the ACC 
field increases. 

Our experimental data are quite similar to the patterns observed on the seminal obser- 
vations of vortices in Bose-condensate gases [s], ^. The clear difference from our work and 
previously related work is that our vortices observed after TOF do not seem to be symmetri- 
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FIG. 6: Mean number of vortex observed in condensed clouds as a function of excitation amplitude 
range for a fixed excitation period of 20ms. After excitation atoms are left trapped for additional 
20ms. 

cally distributed along the condensate j^. Rather, from shot-to-shot they appear at different 
relative spatial positions. In the case of a single vortex, we have observed its position at 
the center of the condensate or towards the edges of it. A possible reason for these unusual 
patterns observed is that we wait only 20ms with the atoms held in the magnetic trap prior 
to release and the well-known regular structures typically need more time (on the order of 
hundreds of milliseconds) to be formed. In other hand, our results resemble the results ob- 
served in Ref . [l^ . Possibly, the similarities can be regarded to the stochastic nature of the 
excitation. In their case, it is due to the relative phase of the merging condensates. In our 
case, two possibilities arise. First, it is possibly due to the fact that one cannot say in which 
direction of the oscillatory motion of the condensate the vortex is excited and it is likely 
to vary from shot-to-shot. Second, if we are really imprinting Berry phase on the conden- 
sate, the decay mechanisms might be responsible for these observations and unusual spatial 
distributions. Nevertheless, to imprint Berry phase on the BEC, the protocols described in 



16 



17| propose a continuous exchange from quadrupolar to hexapolar potentials so break- 
ing the time-reversal symmetry of the processes of winding and unwinding the condensate. 
That procedure gives the sample net angular momentum and promotes the appearance of 
vortices. In addition, in both vortex pump protocols, the inversion of the magnetic bias field 
is an essential requirement. The non-desired effects due to that change like spin- flips and 
subsequent loss of atoms from the trap are overcome by adding an extra optical potential. In 
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our case, there are no such changes. The potential remains harmonic, without inversion of 
the bias field as it is shown in the simulation of Figl2l Increasing the time between the end 
of excitation and releasing the atoms from the trap would help to discriminate the dynamics 
of the vortex patterns observed. In particular, the observation of heating might indicate 
the decay mechanism of the vortices observed. For this set of measurements, we did not 
investigate for much longer times of evolution of the BEC, but for slightly longer times (~30 
ms) no considerable heating is observed. 

In order to try to discriminate which is the possible mechanism of vortex nucleation in 
our experiment we performed 2D numerical simulations. The preliminary simulations con- 
sisted in considering the time dependent Gross-Pitaevskii equation (GPE) plus an external 
potential similar to the one generated by our coils. We find a rather complex behavior of 
the condensate dynamics compared to the cases studied so far of vortex nucleation in the 
presence of a weak elliptical deformation of a rotating trap, where vortices emerge in a rather 
straightforward way 22, l23|. This is due to the elongated geometry of our condensate, that 
is subjected to large deformations for our experimental procedure, making the numerical 
analysis of possible formation of vortices or turbulent regimes 2J] more difficult. In our 
case, the effect of the thermal cloud is included in an effective way by means of noise on 
the initial wavefunction. However, we are not able to simulate explicitly the friction be- 
tween the condensed and thermal clouds. In fact, we believe that a possible mechanism of 
formation of the vortices might be related to the relative movement of the condensed and 
thermal components subjected to the external field. This phenomenon is well-known in the 
scope of fiuid interface theory and experiments and is named Kelvin-Helmholtz instabilities 



25(1 . These instabilities occur in the surface of two fiuids that have a relative velocity and 



give rise to the formation of vortices in the interface of them. This phenomenon has been 
studied in mixtures of superfiuid-normal fiuid liquid Helium 26|] and has been theoretically 
investigated in this context [s^]. Nevertheless, the Kelvin-Helmholtz instabilities have never 
been directly observed as interface excitations in quantum atomic fiuids, though evidence 



of Kelvin modes of a vortex line have been observed [28| in the decay of a counter rotating 
mode in the presence of a quantized vortex in a Bose-Einstein condensate. 

We claim that this is a possible mechanism of formation of our vortices because, if one 
changes the contrast of our images to highlight the interface between the condensed com- 
ponent to the normal component of the gas, round structures can be seen all around the 
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FIG. 7: Atomic optical density image of the atomic cloud from Figl5]^c) with a different contrast, 
showing round structures around the condensed component. 

condensed cloud, as shown in Figd This image is the one shown in FiglS^^c) where only the 
contrast has been modified. This experimental observation performed just after excitation 
reveals large population of vortex-like excitations at the interface and certainly indicates 
towards the mechanism proposed. 

Currently, more thorough simulations are under way to confirm or deny these hypotheses 
as well as trying to apply the standard Kelvin-Helmholtz instabilities theory to quantum 
atomic fluids in order to justify and clarify the true nature of the vortex nucleation observed 
in our sample. The full theoretical treatment is out of the scope of this communication and 
will be addressed in a future publication. 

An important peculiarity of our experiment is that our way to produce vortices is quite 
different from the conventional introduction of rotation by a "spoon" , which always results 
in a single sense of rotation and therefore the formation of vortices with the same circulation 
sign. For our case, the oscillation in the excitation represents equal chances to form vortex 
(positive-circulation) and anti- vortex (negative circulation). In fact, we believe that both 
types are being simultaneously produced in the experiment. This fact may well be part of 
the reasons for observations of such large fluctuations in the response of the system under 
equivalent experimental conditions. 
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IV. CONCLUSIONS 



To summarize, we produced vortices in a BEC of ^'^Rb atoms by superimposing an os- 
cillatory field to the magnetic trap. By changing the amplitude of the external excitation 
we can observe a regime without vortices and another regime where several vortices are 
nucleated. The mechanism responsible for the vortex nucleation due to oscillations seems 
to be related with the presence of a thermal component and its relative motion with the 
condensed fraction. At present, we are performing numerical simulations for proving this 
hypothesis. 

Another important characteristic of our experiment is that the oscillatory motion condi- 
tions can generate both circulation sign vortices and therefore creating interesting conditions 
for vortex/anti- vortex experiments. A systematic investigation increasing the in-trap evo- 
lution time of the sample should reveal details on the relaxation processes, specially the 
possible vortex-antivortex annihilation that would result in heating or the appearance of 
more regular patterns. 

Finally, as previously shown, increasing the amplitude of the excitation, the number of 
vortices nucleated increases dramatically, opening up the possibility of generating turbulent 
regimes in the quantum fluid. Quantum turbulence was discovered in superfiuid "^He more 
than 50 years ago |29|. However, it has never been observed in a quantum degenerate 



gas. A turbulent regime in such a system would represent a new scenario for studying this 
phenomenon, bringing new insights on it. 
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